Vortex Plastic Motion in Twinned Superconductors 
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We present simulations, without electrodynamical assumptions, of B(x,y,H(t)), M(H(t)), and 
J c {H(t)), in hard superconductors, for a variety of twin-boundary pinning potential parameters, 
and for a range of values of the density and strength of the pinning sites. We numerically solve the 
overdamped equations of motion of up to 10 4 flux-gradient-driven vortices which can be temporarily 
trapped at ~ 10 6 pinning centers. These simulations relate macroscopic measurements (e.g., M(H), 
"flame" shaped B(x, y) profiles) with the underlying microscopic pinning landscape and the plastic 
dynamics of individual vortices. 

PACS numbers: 74.60. Ge 



One of the major unsolved problems in high-T c super- 
conductivity is flux pinning . Not only is flux pinning 
of intrinsic interest, understanding it is necessary in or- 
der to develop technological applications. Since vortices 
are pinned by interacting with defects, a great deal of 
work has focused on naturally occurring defects in high- 
Tc superconductors In particular, it has been found 
that twin boundaries (TBs) [2-19] form easily in the 
high-T c compound YBa, 2 C\x z 7 - x (YBCO). Thus TBs 
have been intensively studied by numerous experimental 
techniques, including Bitter decoration |^], torque mag- 
netometryB magnetization and magneto-optical 
imaging |10[-14|. In spite of this intense and multi- 
pronged effort, there is still controversy concerning the ef- 
fect of TBs (see, e.g., the conflicting statements of ||7|,^5|). 
In particular, the magneto-optical measurements by Du- 
ran et a l JTot and Vlasko-Vlasov et al. [jllj seem contra- 
dictory: |1C|] finds that TBs act as easy-flow channels for 
longitudinal vortex motion, while [TT| finds them to be 
barriers. Moreover, resistivity measurements by Kwok 
et al. suggest that pinning is enhanced within the 
TB. Recent results confirm that the effect of TBs 
on vortex motion has a nontrivial dependence on several 
variables. First, there is an angular dependence, i.e., a 
twin may act either as an easy flow channel or a barrier 
depending on the angle between the flux front and the 
TB. Second, experiments |6|,[?],[l4| indicate that the effec- 
tive pinning of TBs is temperature dependent. A com- 
plete picture is still lacking due to the limitations of the 
techniques used so far (e.g., magneto-optics imaging only 
probes the regime of low fields, while resistivity explores 
only the high temperature regime). 

An alternative technique for investigating the micro- 
scopic behavior of flux in a hard superconductor is to 
use computer simulations (see, e.g., [17-1^], and refer- 



ences therein). Here, we present simulations of the spatio- 
temporal evolution of rigid flux-gradient-driven flux lines 
in a twinned superconductor such as YBCO. We first de- 
scribe a model for vortex-vortex, vortex-pin, and vortex- 
twin interactions. We then investigate how these flux- 
gradient-driven flux lines travel through the sample in the 
presence of both correlated (TB) and uncorrelated (point 



defect) pinning. From the microscopic dynamics of indi- 
vidual vortices, we compute several macroscopic measur- 
able quantities (e.g., B(x,y,H(t)), M(H(t)), J c (H(t))) 
and relate them with pinning parameters: the angle 
between the TB and the sample edge, and the ratio of 
pinning strengths inside and outside the TB, which is 
expected to vary with temperature T . 

Our simulation geometry is that of an infinite slab of 
superconductor in a magnetic field applied parallel to 
both the slab surface and the TBs. Thus, demagnetiza- 
tion effects are unimportant. We also treat the vortices 
as stiff, so that we need to model only a two-dimensional 
(2D) slice (x-y plane) of the 3D slab. One can think of 
our rigid flux lines as representing the "center of mass" 
positions of real somewhat-flexible vortices, and our pin- 
ning in the bulk as representing the average of the pinning 
along the length of the real vortex. For a flexible vortex, 
the bulk pinning can be comparable to the pinning in 
the TB even for a large sample. Vortices enter the sam- 
ple at points where the local energy — as determined by 
the local pinning and vortex interaction — is low. We cor- 
rectly model the vortex- vortex force interaction by using 
the modified Bessel function K\(r/X) with cutoffs at 0.1 A 
and 6A (A is the penetration depth). Thus, the force (per 
unit length) on vortex i due to other vortices (ignoring 



cutoffs) is f™ = Ef=i Jo Ki(\Ti -r s \/\) fy . Here, the 
Tj are the positions of the N v vortices within a radius 



6A, 



(r, - rj)/\ri - r,-|, and / 
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Btt 2 A 3 [0. 

Forces are measured in units of /o, lengths in units of 
A, and fields in units of <E>o/A 2 . Notice that our simula- 
tion correctly models the driving force as a result of local 
interactions; no externally imposed "uniform" current is 
applied to the vortices. Further details of the simulation 
are presented in p9| . 

Uncorrelated Pinning. — The most important imper- 
fections in YBCO single crystals are random microscopic 
defects such as oxygen vacancies in the Cu02 planes. In 
order to model such a large number of defects, we divide 
our system into a very fine (e.g., 1000 x 1000) grid, where 
each grid site represents a pinning site. The density of 
pinning sites, n p , is between 150/ A 2 and 496/ A 2 for all the 
runs presented here. These densities are within experi- 
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mentally determined values |q). The maximum pinning 
force ft~E on each grid site (/, m) is randomly chosen with 
a uniform probability distribution, from [0, f p ], where 
f p is the maximum possible pinning threshold force (or 
strength). If the site is in a TB, then the maximum pos- 
sible pinning threshold force is f p B . When the net force 
acting on a vortex i located on grid site (I, m) is less than 
the threshold force then the vortex remains pinned 
in that grid site. However, if the force acting on that vor- 
tex becomes larger than the threshold force ft^, then the 
effective pinning force drops to zero and the vortex moves 
continuously until it encounters a pinning site that has a 
threshold force greater than the net force acting on the 
vortex. Thus the pinning force on each grid site acts like 
a stick-slip friction force f" p = - ff ct 0(/#£ - /f ot ), 
where ff ot is the net force on the vortex (due to other 
vortices and possibly a TB) and 6 is the Heaviside step 
function. In other words, f- p = - ff ot if /f et < 

and = if /f ct > ftp 

Twin Boundary Pinning. — We have tested a very large 
variety of models for TBs, and compared simulation re- 
sults to experiment. For the sake of concreteness, we 
focus here on the model which was most consistent with 
experimental data; the results of the other models will be 
discussed elsewhere. In this model, a TB consists of an 
attractive well containing pinning sites that have a dif- 
ferent maximum threshold force (or strength) /™ than 
that for sites outside the TB, f p . The ratio f p B /f P is 
expected to vary as a function of temperature T because 
of the reduced dimensionality of the thermal fluctuations 
for vortices in a TB The ratio is small (i.e., pinning 
in the twin is weak) for low T, and it is larger at higher 
T (see, e.g., page 1321 of ||). By varying this ratio, we 
can mimic some of the effects due to temperature. 

For all the results presented here, wc model the 
attractive TB well as a parabolic channel with a 
width 2£ TB . The attractive force on the ith vor- 
tex due to the well of the £:th twin is thus f^ TB = 
f TB (dl B /e B ) - dJ k B )/X) % , where dJ B is the 

perpendicular distance between the ith vortex and the 
fcth TB, and / TB is the maximum force that the well 
of the twin exerts on a vortex (i.e., the force needed to 
escape the TB). The pinning inside the twin channel is 
modeled identically to the pinning in the rest of the sam- 
ple, except that it has a different maximum strength, /™ 
instead of f p . This microscopic model of TB pinning is 
very similar to the one inferred from recent measurements 
in M , where it is found that the TB "channel" has strong 
"depth" variations. 

Dynamics. — The overall equation for the overdamped 
motion of a vortex moving with a velocity v and sub- 
ject to vortex-vortex and vortex-pinning forces (due to 
intrinsic pinning f^ p and TB pinning /^ TB ) is U — 

fvv + f vp + f vTB = 7?Vj ) where f . = (f „„ + fvTB) [! _ 

" ff B ~ fDl and we take r, = 1. 



Many parameters can be varied, making the system- 
atic study of this problem quite complex. Here we chose 
to vary only the two most critical variables, the angle 
between the TB and the flux flow, as stressed in [^2[-pT[ , 
and the ratio f p B /f P - We find the critical current, J c , 
for an untwinned sample with f p = 0.2/o and A = 1400A, 
to be 1.26 x 10 6 A/cm 2 , which is reasonable for YBCO 
at low T. Thus we varied f p between 0.15/o — 0.4/q. 
A more detailed investigation with different TB widths, 
TB depths, intrinsic pinning strengths and different TB 
potentials will be presented elsewhere. For the runs pre- 
sented here, we fixed the depth of the TB channel by 
setting the maximum TB force f TB to 0.6/o, and the 
width of the channel 2£ TB to be approximately 0.7A. 

Sample. — We begin by studying in detail how vortices 
entering a sample interact with a single TB. The actual 
sample region is heavily pinned, and covers the region 
with (x,y) coordinates < x < 66A and < y < 55A 
(Fig. 1). The number of pinning grid sites in the sample 
is I x 10 6 = 833, 333. The linear size of each pinning grid 
site is d p = 66A/1000 = 0.066A, and the pinning density 
is n p = 1/dp = 230/A 2 which, for a computer simulation, 
is extremely large. Fig. 4 shows results for other samples. 

Flux Density Profiles and Current Flow. — Figure 1 
shows the computed flux distribution for f p = 0.15/o 
with fJ B /f p = 3/4 (a,b,c) and fJ B /f p = 1/4 (d). The 
TB is at an angle 8 = 60° measured from the horizon- 
tal sample edge located at the bottom, where flux en- 
ters (note the similarity to the experimental situation in 
fj"Tf). Fig. 1(a) shows the positions of the vortices after 
the external field has been ramped up from zero. Two 
important features stand out. First, we see that the vor- 
tices "pile up" on the side of the twin where the flux 
enters. Second, on the other side of the TB we see the 
"shadow" effect; the density of vortices is lowered. These 
features of the vortex plastic flow are more explicit in (b) , 
which shows the induction B along the horizontal line in 
(a). When the external field is ramped down to zero we 
also observe some shadowing and pile up effects, but they 
are weaker since the number of vortices in the sample is 
now smaller. Both of these effects for the ramp up and 
down cases are clearly observed in the experiments of 
pT| . Since we know the exact position of the vortices, we 
can calculate the B field due to the vortices, and hence 
the current flow J{x, y, H(t)) (contours of constant B) 
near the TB. Fig. 1(c) shows the current flow for the TB 
in (a), while (d), corresponding to lower T, presents a 
"flame" -shaped current flow similar to the ones observed 
in several experiments jio] . 

In Fig. 2 we show the positions of the vortices after 
the field has been ramped up from zero for several cases. 
In the top panels, the ratio f p B /f P is held fixed, and the 
angle 8 is varied; we see that the 8 = 20° TB is a barrier 
to flux flow, while the other angles are easy-flow chan- 
nels. Our model, like experiments (e.g., |p^| , |l4]| ), exhibits 
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TB vortex flow which has an angular dependence. In 
the lower panels of Fig. 2 the ratio /™// P is increased. 
We still observe an angular dependence, but the angle at 
which a twin becomes an easy-flow channel has increased. 
As the ratio /™// P increases, the angle a between the 
TB and the direction of vortex flow must decrease if the 
twin is to be an easy flow channel. 

Figure 3 summarizes results obtained from many runs. 
Fig. 3(a) presents the "vortex plastic motion" phase dia- 
gram showing the region where the TB is an easy-flow 
channel (barrier) for longitudinal vortex motion, indi- 
cated by open circles (solid squares). The boundary, in- 
dicated by grey diamonds, corresponds to slow vortex 
motion. This phase diagram quantitatively summarizes 
how 9 and fp B /f p affect the vortex motion. As pinning 
in the twin increases, we expect the flux- front velocity 
to decrease. This observation is quantified in Fig. 3(b). 
There, v TB (v) is the velocity of the flux front inside (out- 
side) the TB, measured along the direction perpendicular 
to the sample edge. Figure 3(b) shows v TB jv decreasing 
monotonically for increasing fJ B /f P - Monitoring the dy- 
namics of individual vortices, we observe that the longi- 
tudinal vortex transport in the TB is plastic (many vor- 
tices "tear away" from the flux lattice when v TB jv ^ 1) 
and proceeds via two mechanisms: (i) the actual mo- 
tion of vortices along the TB, and (ii) vortices from the 
bulk (and close to the TB vortex front) falling into the 
TB. Finally, our phase diagram [Fig. 3(a)] suggests that 
the contradictory experimental results in pO| , pl]| can be 
understood as arising from probing different parameter 
regimes. 

Refs. [p^H made the interesting observation that flux 
enters the sample in a "flame" pattern around the TB. 
jl4j suggests, using current-conservation arguments, that 
this happens whenever the twin acts as easy-flow channel. 
We also observe flame patterns (see Fig. 4(a)) whenever 
the twin is an easy-flow channel. The flame is caused by 
vortices escaping from the TB, not gliding along it. 

Magnetization. — One clear advantage of our simula- 
tion is that we can obtain direct spatio-temporal infor- 
mation on the distribution of individual flux lines inside 
the sample. This is quite difficult to obtain experimen- 
tally, especially for bulk samples. What is typically mea- 
sured in experiments is the average magnetization M over 
the sample volume. In our simulation, we thus calculate 
M = (AirV)- 1 j(H - B) dV . In Fig. 4(b) we show par- 
tial magnetization loops for a sample where the TBs are 
perpendicular to the flux flow (d), one where they are 
parallel to the flux flow (e), and a sample without any 
twins at all. Note that Fig. 4(d,e) are two particular (and 
extreme) cases considered in the phase diagram Fig. 3(a): 
9 = 0° and 90°. For these runs the sample is 42 A x 42 A, 
fJ B /f P = 0.1, and f p = 0.2/q. 

The loop where the twins are perpendicular to the flux 
flow is clearly wider than the other two. This is because 
a twin is a barrier to perpendicular vortex flow since the 



vortices would have to cross the TB, which is a parabolic 
channel. The loop where the TBs are parallel to the flux 
flow is narrower because these TBs are easy flow chan- 
nels. Also, we repeated the above experiments at a higher 
ratio, fp B /f p — 1.5, and found that the perpendicular 
TBs arc still barriers to flux flow, i.e., the magnetization 
is still higher than the untwinned case. The parallel twins 
no longer act as easy-flow channels and thus the magneti- 
zation is nearly identical to the untwinned case. All these 
results are qualitatively similar to those in Ref. [Q; the 
minor differences are most likely due to finite size effects. 

To summarize, we have performed simulations of flux- 
gradient-driven vortices in twinned superconductors with 
a very large number of (~ 10 6 ) point defects. Experi- 
ments can be explained by using a simple model for a 
twin boundary: a parabolic well containing pinning sites 
which are different in strength from the rest of the sam- 
ple. In agreement with experiments, we show that the 
twin may act as either a barrier or an easy flow chan- 
nel (producing "flame" -type flux patterns), depending on 
the angle between the flux flow and the pinning ratio, 
fp B /f P - We quantify these results in detail and pro- 
vide a phase diagram for the plastic motion of vortices 
[Fig. 3(a)]. When the twin is a barrier, it acts to guide 
flux motion resulting in the pile up and shadowing effects. 
We also compute magnetization loops for samples where 
the twins are parallel and perpendicular to the flux flow 
and obtain results similar to recent experiments 
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x (in units of X) x 
FIG. 1. (a) Snapshot during the initial ramp-up phase of 
the top view of the region where flux lines, indicated by dots, 
move, (b) B(x) along a 1. OA-wide horizontal cut through 
the sample centered on the horizontal line y = 36A indicated 
in (a). We calculate B(x,y) at each point due to all the 
vortices and then find B(x) by averaging over a strip, i.e., 
B(x,y = y ) = (Ay)" 1 J^ 0+Ay dy B(x, y). This B{x) clearly 
shows the "shadowing" and "pile up" effects discussed in jy] . 
(c,d) Current flow (contours of constant B) around the TB in 
(a) for several pinning strengths. In (c) f p /f p = 3/4, as in 
(a), more effective source of pinning, and the TB is a barrier 
to longitudinal vortex motion. In (d), f p B /f P = 1/4 and the 
TB is an easy-flow vortex channel. 



8 = 20, f™/f„ = 1/4 


6 = 45, fJ7f p = 1/4 


8 = 85, f™/f„ = 1/4 






= 20 f™/f, = 1/2 


6 = 45. I 'll = 1/2 


8 = 85, f™/f„= 1/2 



FIG. 2. Top view of the region where flux lines, indi- 
cated by dots, move in the ramp-up phase. We have varied 

8 and fJ B /f P , and here we present some of our results for 

9 = 20° (left), 45° (center), 85° (right), and f p B /f p = 1/4 
(top), 1/2 (bottom). For the top panels the TB acts as an 
easy-flow channel, except = 20°, while for the bottom ones 
only the sample with 8 — 85° exhibits easy-flow channel be- 
havior. Here f p — 0.15/o and n p — 230/A 2 . 
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FIG. 3. (a) Phase diagram /™ / f p versus 6 showing the 
region where the TB is an easy-flow channel (barrier) for 
longitudinal vortex motion, indicated by open circles (solid 
squares). The boundary, indicated by grey diamonds, cor- 
responds to slow vortex motion. The ratio v TB /v is ~ 1 
(black squares), between 1 and 1.25 (grey diamonds), and 
larger than 1.25 (open circles), (b) v TB /v versus f p B /f P for 
= 35°, 70°, 80° (upper curve). 
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(e) " (d) 

(a) B z {x,y) along a TB at 6 = 85° with 
f P B /f P = 0.05, /„ = 0.4/o, and n p = 152/A 2 . B z (x,y) ex- 
hibits a "flame" front profile. The same-color contours of 
constant B z indicate the current flow around the twin. The 
images 4(a,d,e) show B z (x,y), which is obtained by approx- 
imating the field due to each vortex by a monopole located 
1A below the surface and computing B z (x,y) at z = 0.6A 
above the surface (see Ref. pQ]). (b) M(H(t)) for heavily 
twinned samples; the TBs are 5.6A apart. TBs, indicated by 
black lines in (d) and (e), act as barriers (for TBs perpendic- 
ular to the vortex flow as in (d)) or as easy-flow channels (for 
TBs parallel to the vortex flow as in (e)). For comparison, 
M(H(t)) for an untwinned sample is also shown in (b). In 
(b-e), ff B /f p = 0.1, f p = 0.2/o, and n p = 496/A 2 . (c) B{x) 
for the sample shown in (d). Notice the sudden decrease in 
B(x) (leading to high-currents) at each of the TBs. This "ter- 
raced" B(x) profile produces a periodic array of alternating 
high and low current regions (see Ref. (il)). In (c,d) the TBs 
act as a series of "flux dams" . 
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